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Abstract 

We study several models of random geometric subdivisions arising 
from the model of Diaconis and Miclo (2011). In particular, we show 
that the limiting shape of an indefinite subdivision of a quadrilateral is 
a.s. a parallelogram. We also show that the geometric subdivisions of a 
triangle by angle bisectors converge (only weakly) to a non-atomic dis- 
tribution, and that the geometric subdivisions of a triangle by choosing 
random points on its sides converges to a "flat" triangle, similarly to 
the result of Diaconis and Miclo (2011). 

Keywords: barycentric subdivision, geometric probability, Markov chain, it- 
erated random functions. 
AMS subject classification: 60J05, 60D05. 

1 Motivation 

The aim of this paper is to consider several models involving random subdivi- 
sion of geometrical objects. Markov chains involving geometry and polygons 
have been studied quite widely in the literature, see e.g. [3] and references 
therein; however, the main motivation of this paper comes from the paper by 
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Diaconis and Miclo [2] , who considered the following model, earlier also stud- 
ied by Hough [5]. A non-degenerate triangle is divided by its three medians 
into 6 smaller triangles, and one of these new triangles is chosen with equal 
probability, and (optionally) re-scaled, thus becoming the "new" triangle. 
This procedure is repeated indefinitely It was shown that in some sense the 
limiting triangle will be "flat", that is the largest angle will converge to tt. 

In the current paper we consider several generalizations of the above 
model. In Section [2] we consider subdivisions of a quadrilateral by the lines 
connecting the middle points on the opposite sides; in Section [3] we consider 
subdivision of a triangle using angle bisectors. Finally, in Section H] we con- 
sider a sequence of triangles obtained by randomly choosing a point on each 
of the sides and letting them be vertices of the new triangle. 



2 Random subdivision of quadrilateral 




We are given a convex quadrilateral ABCD. Let E, F, G, H be the middle 
points of the sides AB, BC, CD, and DA respectively. Let M be the point 
of intersection of segments EG and FH. Now we replace ABCD by one 
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of the following four quadrilaterals AEMH, EBFM, MFCG and HMGD 
with equal probabilities. Suppose that we repeat this procedure indefinitely. 
What is the limiting shape of the quadrilateral obtained as the limit of this 
procedure? 

Theorem 1 Under the procedure described above, the limiting shape of the 
quadrilateral will be a parallelogram, and the rate of convergence is geometric. 
(Note that the shape of parallelogram is "invariant" for the procedure). 

Proof. Observe that HF = \{AB + DC), and HM = MF = \HF where 
x denotes the vector x. Also suppose that when we replace the original 
quadrilateral by one of the four smaller ones, we rescale the smaller one 

— * 

twice thus making it bigger; this will not affect the shape. Let u = AB and 

— * 

v = DC be the vectors corresponding to the "horizontal" sides of ABCD, 
and Mi, v\ be the corresponding vectors of the new quadrilateral obtained by 
subdivision. The crucial observation is that 




Let us place all vectors u n ,v n at the origin O and let U n and V n be the 
corresponding endpoints of these vectors. Let X n = U n if U n £ {C/„_i, V n -i} 
and X n = V n if V n {U n -±, V n -i} (exactly one of these two statements must 
be true). Then we see that all points X n lie on the segment U V . Moreover, 
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Xi lies exactly in the middle of U V , X 2 with equal probabilities splits U Xi 
or V X 1 in the middle, etc. As a result, we see that X n a.s. converges to 
a point Xoo which is uniformly distributed on the segment UqVq. Taking 
into account that |£/ n 14| = 2~ n |[/ V r o|, we obtain a deterministic speed of 
convergence of u n and v n towards := OX^: 

\UoVo\ { ._ -> \ \-> -> n ^ l^o^ol 

< max{|M n - Moo|, \v n - Mod} < — ^— ■ 

The analogous statement holds also for the "vertical" sides corresponding to 
AD and BC, thus yielding the desired convergence towards a parallelogram. 
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Unlike the median-subdivision model considered in [2J, suppose that we 
subdivide the triangle ABC by the three angle bisectors which intersect the 
sides AB, BC, CA at points E, F, G respectively, and let M be the point 
of intersections of all angle bisectors, the centre of the inscribed circle. The 
replacement procedure now states that the triangle ABC is replaced with 
probabilities 1/6 by one of the following triangles: ADM, DBM, BEM, 
ECM, CFM, FAM. As before, the object of interest is the shape of the 
limiting triangle obtained by indefinite repetition of the above replacement 
procedure. 

It turns out to be convenient to work with the angles of the triangle. If 
the original triangle has the angles (a,/3, 7), a + B + 7 = n, then the new 
triangle will have one of the following six sets of angles: 
a 8 a + 8\ fa „ 7 a + 7 



2' ' 2' 2 J ' \2 ' 2' 2 
^7 P+j\ fP , a a + P\ ,,, 



2 U 2' 2 y ' ^2'" ' 2' 2 
Obviously, we cannot expect convergence almost surely for this procedure 
(since e.g. any angle can be halved on the next step with probability 1/3). 

Observe that we can generate the sequence of triangles by always choosing 
the left-bottom triangle, that is by using the mapping (a, /3, 7) — >■ {%, 7 + f , ^2 ) ' 
and then performing a random permutation of the set of three newly obtained 
angles. Formally, let (cx n ,j3 n , j n ) denote the set of the angles of the n-th tri- 
angle, then 

in \ ( ^n Pn &n Pn 
(a n +l,Pn+l>7n+lJ = &n — , 7„ + — , 



2 2 2 

Pn ^n Pn\ / <->\ 

■y^-^-T'— 2- ) (3 ' 2) 

where a n = {p^, <j( 2 \ a^ 3 \ a^ 5 \ a*- 6 )} is a random permutation of the 
set of three elements; cr n ((a, b, c)) takes one of the six possible values 

{(a,b,c), (b,c,a), (c,a,6), (c,6,a), (6,a,c), (a,c,6)} 
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with equal probability, and o"„'s are i.i.d. The advantage of this method is 
that the distribution of the set of the angles of the n-th triangle is completely 
symmetric with respect to exchanges of its components, even though the 
components are not independent. Unlike the case of barycentric subdivision 
studied in [2J, only the weak limit exists in our case. 

Theorem 2 The sequence (a n , f3 n , j n ) converges in distribution to some limit. 

Proof. Let S = {(a, (3, 7) : a, /3,7>0, a + (3 + 7 = n} be the 2-simplex with 
the standard Euclidean distance. Let fi(u), u 6 S, i = 1, . . . , 6 be the set of 
functions given by (13. ip . Let u = (a, (3, 7) and v — (a + x,fl + y,j + z), so 
that x + y + z = 0. Observe that all /j are Lipschitz and that 



Ej \fj( u ) - fj( v )\ _ 1q / V 1 + V 2 + [v + A 2 + xz \ i x2 + z 2 + [y + z] 2 + xy] 
^ 111. 91 1 



1=1 



\u — v\ y 4(x 2 + y 2 + z 



2\2 



1 (y 2 + 3^ 2 + 3yz) (3y 2 + ^ 2 + 3yz) 
~ 2 ° g 16(2/ 2 + <2 2 + l/^) 2 

1 (3y 2 + 3^ 2 + 3^)(3y 2 + 3^ 2 + 3^) . 3 
^2 l0g 16(, 2 + , 2 + ^) 2 = 1 ° g 4 <0 - 

The identical bound holds for % = 3,4 and z = 5, 6. Therefore, if Z u denotes 
a set of the angles obtained from u by random subdivision, we have 



E 



log 



\Z U 


~~ Z v \ 


\u 


— v\ 



E i log m^M< log f <0 . 

^-^ o \u — v\ 2 

i=i 1 1 



Therefore, by Theorem 1 from [TJ, see also the proof of Lemma 5.1 in [2J, the 
mapping u — > Z u is ergodic, that is there is a (unique) probability measure 
v on S such that for any starting configuration u = (a,(3,j) G S we have 
— > ^ weakly. (Here stands for the superposition of n i.i.d. mappings 
Z.) m 

Now let us try to get a handle on the distribution of the limiting triangle. 
For the purpose of simplicity, and without loss of generality, assume from 
now on that a n + (3 n + 7„ = 1 (as opposed to n). 
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We already know that (a n , (3 n ) converges to some pair (a, (3) in distribu- 
tion. Since a is bounded and f3 and 7 := 1 — a — (3 have the same distribution 
as a, and a (3, {3j, ■ya all have the same distribution as well (from the sym- 
metry), we conclude 



E [a + + 7] = 1 
E [a + p + 7] 2 = 1 



Ea = E(3 = E7 = -, 
3 

6E [a/3] = 1 



3E« 2 



Also from (13.21) we have 



Ea 2 



This in turn yields 



Ea +E 



7 



(3^ 



1 



E(a + ^) 2 



Ea 2 = J, E[ap] = —, Var (a) = — , CWa, /3) = - — 
7 L J 21' v ; 63' v ' ' 63 

which sheds some light on the distribution of a and the dependence between 
a and (3. 

A more interesting and subtle statement about the joint distribution of 
a, P, 7 is the following 

Theorem 3 Let c±, C2, C3 6e some real numbers, not all of which are 0. Then 
distribution of the random variable C\a + c 2 /3 + c 3 7 c?oes not have atoms. 

In fact, we conjecture that the distribution of v is continuous on the 
simplex, and so are the marginal distributions, e.g. the distribution of a. 
Numerical simulations suggest that the pdf of a randomly chosen angle of 
a limiting triangle looks like the one shown on Figure [TJ which is obviously 
quite non-trivial. 

Due to the symmetry of the triple (a, (3, 7) with respect to permutations 
and the fact that a + (3 + 7 — 1 > Theorem [3] follows immediately from the 
next statement. 



Lemma 1 For any c G K and 16I, P(cc + cp = x) = 0. 
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O 1 2 3 

x 



Figure 1: Density of a randomly chosen angle. 



Proof. Equation ( 13 .2p together with the weak convergence of (a n , (3 n , 7 n ) yield 



[a + c/3 = x) = -<¥[- + c 
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a 



1 - & 
1 — a 

2 



.... a a + 8 

+ P - + c - 

1 2 2 



1 - & 



+ P 
+ P I c 



. .... a a + /3 
' 1 r F ' C 2 + 



a + /? 
c — - — = x 



a + 8 
H — = x 



(3.3) 



In particular, 



(a = x) = - {P(a/2 = x) + P(l - a - 8/2 = x) + P([a + /3]/2 = x)} 



(3.4) 



and if we set x = 1, then 



(a = 1) = - {P(a/2 = 1) + P(l - a - 8/2 = 1) + 



7]/2 = l)} 



3 1 



(a = J = 0) + 0} = ip(7 = 1) = ^P(a = 1), 



yielding 



P(a = 1) = 0. 



(3.5) 
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Figure 2: Sample density of v on simplex S. 
Rewriting (13. 3p for c ^ { — 1, 1/2, 2} we have 

P(a + cP = x) = l - |p (a + -^—fi =?J + P (a + = 

+ pfa + -^^=?^ +p( < a + -J— ) (3.6) 



1 + c / V 1 + c 

where the question marks stand for some non-random numbers. We can 
make sense of the expression above also for c G { — 1, 1/2, 2} if we replace the 
fractions equal to infinity by zeros, due to the symmetry between at and j3. 

For each c, let x*(c) be such that P(a+c/3 = x*(c)) = p*(c) = max^R P(a+ 
c(3 = x); clearly such x*(c) must exist. Suppose that the statement of the 
lemma does not hold; then p* = sup ceK p*(c) > 0. From (13.61) and the sym- 
metry between a and ft it follows that 

p*(c) =F(a + c(3 = x*(c)) <\{p* +p* +p* + p*(c + 1) + p* + p*} 

6 

= i{5p*+p*(c+l)}, 
6 
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hence 



p*(c + 1) = p*(c) — 6e whenever p*(c) > p* — e. 



Fix a very small e' > and let c be such that p*(c ) > p* — e'. Then for 
<k — c + i, i — 1, 2, . . . , we have > p* — 6 V. 

Let A(c) = : a + c/3 = x*(c)} and Z N = l A(co) + l A(ci) H hlA( Cjv _i)- 

On one hand, 



N-l 



EZ N = J2 F ( A (d)) > N P* ~ 0-2 • Q N e' 



i=0 



since P(A(q)) = p*(ci) > p* - 6V. On the other hand, for 1 < n < iV, 



JV 



E Z^v = F ( Z n >i)<n-l+ F(Z N > n), 



thus 



Suppose that 



¥(Z N >n)> Np* -(n-l)- 0.2 • 6 N e'. 



Np* > (n - 1) + 0.2 • 6 V. 



(3.7) 



Then there is a subset of A(c ), . . . , -A(cjv-i) containing n distinct elements, 
say A(cjJ, . . . , A(cj n ) for which the probability of n r r l n=1 A(c im ) is strictly pos- 
itive. We are going to make use of 

Claim 1 Suppose that z\, Z2, ■ ■ ■ ,z n , n > 3, is a collection of n distinct real 
numbers. Let B = f)™ =1 A( Zi ). IfP(B) > then the sets Afa) \B, % = 
1,2, ... ,n, are disjoint and as a result 

(71 \ to n 

\jA( Zl ) = P(B) +J2P(A( Zi ) \ B) = J2P(A( Zi )) -{n- l)P(B). 
7=1 / 7=1 1=1 

Proof of the claim. Since P(B) > 0, the system of equations 



a + zif. 
a + Z2I 



x*(z 2 ) 



k a + z n (3 =x*(z n ) 
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must have a solution in (a, (3). This yields that for all distinct i, j G {1, . . . , n} 
we have 

- = x*( Zi )-x*( Zj ) (3 g) 

Zj 

and thus for any three distinct i,j, k 

X*(Zj) - X*(Zj) _ X*{Zj) -x*(z k ) 
Zi Zj Z{ Z]$ 

which in turn leads to u G A(zi) D -A (2?) =>■ oj G A(zfc) for all k =>■ G S. 
Now the second statement of the claim is trivial. ■ 

From Claim Q] it follows that for B = n^ =1 A(c im ) 

p( m > -i + ELiP(^(0) > p * _ 6V _ J_ 

n — 1 n — 1 

Let e" = 6V + From (ESD we have that if we set y = x *^~ x *^ then 

71-1 M " tf C^-C^ 

P(a = y) = P(/3 = y) > P(£) > - e". At the same time, from ([33]) it 
follows 

- e" < P(a = y) = - {P (a = 2y) + P (a + /3/2 = 2(1 - y)) + P (a + = 2y) } 
= ^{F(a = 2y)+F(a + P/2 = 2(l-y))+F(a = l- 2y) } . 

Since each of the expressions on the RHS does not exceed p*, we have 
min{P(a = 2y),¥(a = 1 - 2y)} > p* - 3e". W.l.o.g. assume that 2y > 1/2. 
Reiterating (E3D and recalling ([33]) (so that P(a > 1) = P(a = 1) = ) we 
obtain 

p* - 3e" < P(a = 2y) 

= ~ |p (a = Ay) + P (a + | = 2(1 - 2y)) + P {a = 1 - 4y) 

< ± {0 + P (a + 0/2 = 2(1 - 2y)) + 0} < y 

leading to a contradiction, provided e" < 2p*/9. To finish the proof, we 
have to demonstrate that there exist n, N, e' instantaneously satisfying this 
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condition as well as (J377J)- Indeed, fix an n so large that n — 1 > 9/p*. Next, 
let N be larger than (n — Finally, set 

, min {Np* — (n — l),p*/9} 
It is easy to see that (13. 7p is fulfilled and moreover 

Ti-l 9 9 9 



4 Random subtriangle of triangle 

jo 



B^L^ 










\b 



Now on each side of the triangle we randomly (independently and uniformly) 
choose points A%, B\, C\. The new triangle is now formed by these three 
points {A\B\C\). Repeat this procedure indefinitely. What is the limit of 
the shape? 

Let a = \BC\, b = \CA\, and c = \AB\. Also let £ a = \BAt\f\BC\, 
£b = | C -Bi | / 1 C^4 1 , and £ c = |ACi|/|AB|. According to our assumption, 
£i, i = a,b, c, are independent uniform [0, 1] random variables. Also let 
a% = \BiC\\, b\ = \C\Ai\, and c% = \AiBi\ be the sides of the new obtained 
triangle. 
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Theorem 4 The limiting triangle shape is a.s. flat (i.e., maximum angle 
converges to ir). Moreover, for any c < 1 + ^jp — ^ = 0.365 . . . 

Vn < e~ cn for all sufficiently large n 

where y n is the ratio of the height of the triangle corresponding to its largest 
side and the length of this side. 

Similarly to [2j, let us rescale the triangle such that the largest side's 
length (say, AB) is 1, and fit this side on the coordinate plane such that 
A = (0,0), B = (1,0). Also suppose that C lies in the upper plane and 
\AC\ > \BC\. Let C = (x,y), then x G [1/2, 1], y > 0, and the pair (x,y) 
completely characterizes the shape of the triangle ABC. We also have 

A x = (1 - (1 - x)£ a , yQ , B 1 = (x(l - y(l - &)) , C x = (£ c , 0) . 

The new side lengths will then satisfy 

a? = [x(l-6)-e c ] 2 + b(l-6)] 2 

c\ = [(1 - (1 - x)Q - x{l - &)] 2 + W ~ia- 6)] 2 

and according to the standard formulas, the area of the triangle with vertex 
coordinates Ai, Bi, and C\ equals 

A := | \tMc + (1 - 0(1 - 6)(1 - Co)] 
Thus the new value of y is now 

2A q fe e e + (i-e t )(i-&)(i-ec) 

max{ai, 6i, Ci} 2 max{a 2 , 6f , cf } 

and also 

I 2max{I 1 7iV if Cl = min{ai,6i.ci}; 

X l = \ 2max{t"ciP if fe l = ^{ai, h, d}; 
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which can be summarized as 



Vi = y-r, 

_ a\ + b\ + c\ - 2 min{a 2 , bj, c\} 
Xl 2m&x{al,bl,cl} 

where 

R 



r = r(x,y) = 



max{a 2 , bf, c[} 
and 

R = Rita, &) = Ubtc + (i - 0(1 - - Q- 

Observe that the denominator of r(x, y) is increasing in y; hence r(x, y) < 
r(x, 0). 

From now on let us assume y = 0. For y = 0, we have that points A±, 
B\, C\ all lie on the horizontal axes with coordinates fi = 1 — (1 — x)£ a , 
v = x{\ — and £ c respectively; we always have < v < ji < 1 since 
v G [0,x] and fi G [x, 1]. Consequently, 

S(x;£a,tb,tc) •■= max{ai,6i,ci} 

A* - if £ c < ^; 
max{/i -£c,V-v,tc-v}={ \i — v if i/ < £ c < \i\ 



i c -v if £ c > /i. 



Therefore, 

R{ia-,ib-,ici 



E[r(*,0) / ^"^) 112 de c = / 1 + / 2 + / 3 

Jo [maxlai^^cilj 2 

r R(Za,Zb,Zc) Ar f 1 gfc^&jcK , f'R^a^b^c) 



Easy algebra gives 

A = (6, + 6,-1) hi 

6* + 1 - a 



/ / 6x + (l-a;)(l-ea) \ 



2 



/, = «. h ( , ': X+ Sf + W1 " 6)(1 " x) 



£ b X + (1 - x)(l - i a ) J 1-X + &X 
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hence 

E [r(x, 0) | 6] = & + 6-1) In (^£±f|) + ~ + * " ~ 

- Z b )(x\z a + & - l) + 6,(1 - 2a) - l) 

(1 L ■ C„r){\ x ■ 

As a result, 

E[r(x,0)\Q = [ E[r(x, 0)|£ a ,&]d& 
Jo 

+ 1 - £,(1 - x)(2x + 3 - 6(2 - x))) 



2^(1-^ + 6^) 

+ (1 - 2&)(1 - 6 + M((l - a: 2 ) log(l - x) + x 2 log(l - 6 + M)] 



and thus 



E[r(a;,?/)]<E[r(x,0)]= / E [r(x, 0) | £ ] d£ a = I. 

Jo 

Let (xo,yo)j ( x i,Vi), ( X 2,V2), • • • be the sequence of coordinates corre- 
sponding to the sequence of subdivided triangles. Assume y > 0, then with 
probability 1 we have y n > for all n. Let T n be the sigma-field generated 
by {(x i: yi), % — 0, 1, . . . , n}. We have just established that 

E [y n +i I F n ] = y n E [r(x n , y n )\ < y n . 

Thus y n G [0, 1] is a supermartingale which must converge a.s. to some limit 
Voo . Let E = {£„ < 0.1, & > 0.9}, then P(£) = 0.01 > 0. On the event E 
we have C\ > 0.8, hence 

r( y) = 6-[66] + (i-6)-[(i-6)(i-6)] < 2_ai < i 

max{ai, &i, ex} 2 cf 3 

Consequently, P(y„+i < y n /3 | J-" n ) > 0.01 which implies y^ — a.s. 
Note that 



E [log / / / log(i2)d£ c d&d£o = — -- = -1.57... 



J0 JO 



9 3 
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and 



E [log S^x; = / d£ a / d& 

Jo Jo 



1 /•! 



log(p - Q di c + / log(/i - i/) d£ c 

I! J V 

5 x 3 log(x) + (1 — x) 3 log(x) 

6 3x(l — x) 

log 4 - 5 . 

> — = -0.602- ■•=: -k 4.9 

6 



+ / io g (e c -v)de 



i-i- 



Observe that 



\ogy n - logy n _i = logr(x n _i,?/ n _i) < logr(x n ,0) 

= log R (e\ C\& n) ) - 2 log S (x n _ i; tfO, &\ £<») 

= p n + 2cr n 

where p n = logi? fe' 1 ', <^ n ' > , ^c™' ) 1 are i.i.d. with mean ^ — | and a n = 

- log Si %n—i'i £,a ^' £c ^ ) some J 7 re adapted random variables respec- 
tively. Since < S'(-) < 1, we have <r ra > 0, and also E (cr n | J- n -i) < K due 
to (14. 9p . Additionally, we have for any positive z 



PK > z | J- n _0 = P (s (x n _ i; £\ d n) ) < e~ 2 | J- n _ 

< P^ - i/ < e~ z | 7"„_i) = P(x n _!6 + (1 - x n _i)(l - £«) < I ^n-i) 

< P(l{ a!B _ 1 >l/ a }Sn-l6 + 1{x„_i<1/2}(1 ~ #n-l)(l - O] < e ~ 2 I ^n-l) 

<p(|<e-| J-n-iJ = 2e~* 

since 1 — £ a has the same uniform (0, 1) distribution as 

The following statement is related to exponential inequalities involving 
martingales; however since we did not find it in the form we needed, we 
present its short proof later. 

Lemma 2 Let T n be an increasing family of a— fields, and a n be T n -adapted 
random variables, possibly unbounded. Suppose that 

E [a n | F n -i\ < k 
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and for some z > and c > we have 

P(|^n| > z | Fn-i] < e~ cz for all z > z. 



Then 



\r-\n 

limsup — — — - < k a.s. 

n. — ton ri 



Now from Lemma |2] and the strong law of large numbers it follows that 
lim sup < hm sup — — h 2 hm sup — — — 

n— too 71 n— >oo 71 n—¥oo 71 

it 2 log 4 
<Ep n + 2n = 1 — = -0.365 . . . 

This established the result of Theorem 0] 

Theorem 5 x n converges to Uniform [1/2,1] distribution. 

Proof. Since y n — > a.s., and x n+ \ is continuous in y n near 0, it will suffice 
to show that, given y n = 0, x n+ \ has £/[l/2, 1] distribution. This will follow, 
in turn, from the following statement: put /i, u, £ c in an increasing order, and 
denote the resulting values < x^ < x^ < x^ < 1, then x '■= ^pjz^rs nas 
U[0, 1] distribution. Indeed, for any z G (0, 1) we have 

F(x <*) = HX < z, £ c < u)+¥( X < z, v < £ c < fj) +P( X < z, & > pi) 
= (I) + (II) + (III). 

We have 

n\ n\ pi/ ( (3.2 CCZ^ ' ZX 

(/)= / d& / d^ a / i { t , U<Z \^=\ { ^::L )Z 

Jo Jo Jo \i-e 6 d-)-& / [ i — 6g(1 _ z) 



(i: I J-.,-) if x < 2; 

if x > 2, 



(//)=/ d£ a / d^ /' If ^ ^ d£ c = z/2, 



Jv \ 1— 5&(I-5J— Sa» J 



<yIU " ) = So ^ Jo ^ I 1 [ ^Mi^lr^ c - j d ^ 



Sz^+z^a: — 3z' ! :e+,z3; — 3zx+x z -r ^ 

6a(l-a0 X ^ Z ' 

if X > Z. 



{l-x) 2 z 2 
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As a result, F(x < z) — (I) + (II) + (III) = z and thus \ has a uniform 
[0, 1] distribution. ■ 

4.1 Proof of Lemma [2] 

Proof. The result follows from the Borel-Cantelli lemma, once we establish 
that for any e > there is a 5 = 5(e) > such that 



P f ^2 °i > ( K + £ ) ra J < e_<5n - ( 4 - 10 ) 
First, by Markov inequality, for any A > 

P ( J^ai > (K + e)n) = P (e x ^ ai > e An(K+e) ) < e ~ Xn(K+£) 



>E 



,1=1 



-An(K+£)jg 



Now to show (14. 10p by induction, it suffices to demonstrate that there are 
A > and 5 > such that E [e Xan | 7"„_i] < e *("+ e )-« for all n. 
Indeed, (see e.g. [I], Lemma 5.7 in Section 1) 



E 



A fc |cr r , 

~~ kJ 



A 



fc /-oo 



< 



kl 

A* 
fc! 



/C2 



z k + 



fc-ll 



|cr n | > 2! I J* n _i) d^ 



fcz^e^ dz 



(Az) fc /A 



therefore, assuming A < min(c, 1), by Dominated convergence theorem we 
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have 



E [e A<J " | T n - X ] = 1 + AE [a n | J^J + ^ 

fc=2 



< 1 + Ak + J^E 



fc=2 

oo 



A fc K 

k x 



k=2 V 7 fc=2 



A fc E [a k n \F n -i] 
k\ 

J 7 n-l 

(Xz) k 



k\ 



< 1 + Xk + 



A 2 



+ X 2 \e~ z -1-z] 



c(c - A) 



(4.11) 



Now by choosing A > sufficiently small, we can make ( 14. lip smaller than 
e \{K+e)-8 > ^ _j_ ^ K _|_ ^ £ _ ^ £ Qr some £ > Q, thus finishing the proof. ■ 
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